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The structure of the dark energy equation of state phase plane holds important information on 
the nature of the physics. We explain the bounds of the freezing and thawing models of scalar field 
dark energy in terms of the tension between the steepness of the potential vs. the Hubble drag. 
Additionally, we extend the phase plane structure to modified gravity theories, examine trajecto- 
ries of models with certain properties, and categorize regions in terms of scalar field hierarchical 
parameters, showing that dark energy is generically not a slow roll phenomenon. 
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I. INTRODUCTION 

The discovery of the acceleration of the cosmic expan- 
sion has thrown physics and astronomy research into a 
ferment of activity, from a search for fundamental the- 
ories to investigation of predictions relating models and 
the cosmological dynamics, to development of astrophys- 
ical surveys yielding improved measurements. 

The acceleration, or more generally the expansion his- 
tory of the scale factor evolution with time, a(t), can be 
equivalently treated by a dark energy pressure to density, 
or equation of state, ratio w(a) [lj. One model, Einstein's 
cosmological constant, predicts w — —1 at all times, but 
generically the dark energy phenomenon has dynamics, 
a time varying w(a). 

It is important to note that the current epoch of accel- 
erated expansion is very different from the early epoch 
of inflation. We know a priori that dark energy does not 
completely dominate the universe now and we do not 
know a priori that dark energy obeys a slow roll approx- 
imation (in fact we will see it is unlikely to). 

In these senses, dark energy is a more challenging 
phenomenon than inflation. We are faced with three 
"Goldilocks" problems: 1) dark energy is dynamically 
important, but not fully dominant, with ~ 75% of the 
total energy density today, 2) the universe is accelerat- 
ing, as from a component with w < —0.8, but the field 
responsible may not be slowly rolling (unless fine tuned) 
as it would if nearly completely potential dominated, and 
3) if it's not the cosmological constant, what happened 
to the cosmological constant? 

To discover whether the physics is the cosmological 
constant, and to distinguish between alternate theories, 
requires measurement of the possible dynamics. 

Caldwell & Linder Q (Paper 1) investigated the dy- 
namics in the phase plane of w-w' , where w' = dw/dlna, 
for canonical scalar field dark energy, or quintessence, 
finding that this reveals important clues to the nature 
of the new physics. In such a plane, the time or scale 
factor variable is a parameter along the paths of dynam- 
ics. They found distinct structure, categorizing fields into 
those that at early times are locked by Hubble friction 
into a cosmological constant like state, and then move 
away from this (thaw) as the dark energy dominates, 
and those that initially roll due to the steepness of the 



potential but later approach the cosmological constant 
(freeze). In this article some of these results are put on a 
firmer footing, examined in greater detail, and extended 
to models beyond canonical scalar fields, including mod- 
ified gravity theories. 



II. GENERAL DYNAMICS 

We begin with a brief, explicit derivation of the key 
dynamics equation, Eq. JBJ. A scalar field (j> possesses 
both potential energy V((f>) and kinetic energy (l/2)</> 2 . 
The Lagrangian density is of the form 

C={l/2)g^d^d"<t>-V{4>), (1) 

for a canonical, minimally coupled scalar field in a metric 
g M „. The equation of motion for the field, the Klein- 
Gordon equation, 



3H<f> + V,A = 0, 



(2) 



where H = a/ a is the Hubble parameter, follows from 
functional variation of the Lagrangian. (Spatial inhomo- 
geneities should be negligible on subhorizon scales, but 
also see WI CI ) The energy- momentum tensor is gener- 
ated through Noether's theorem and one can identify the 
energy density and pressure: 

(l/2)0 2 + l/(0) ; p=(l/2)^-V(^). 



P 



(3) 



It will be convenient to invert these and write the poten- 
tial and kinetic energies in terms of p and w: 

V = (l/2)(l-w)p ; (l/2)^ = (l/2)(l + «;)p. (4) 

Note that w and p are both functions of time, or scale 
factor. 

To obtain an equation for the variation w', take the 
derivative of V with respect to <p, 

V,<t> = V/j> = (1/2)[(1 - w)p - pw]/[(l + w)p] 1 ' 2 . (5) 



Employing the continuity equation p 
H = d\na/dt, one obtains 



-3(l-u; 2 )-(l-w)(l+w) 1/2 v /(3/87r)^(a) 



3Hp(l + w) and 
MpVs 



V 



(6) 

where fl,p(a) = 8ir p / '(3H 2 Mp) is the dimensionless dark 
energy density, and Mp is the Planck mass. 
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A. Distinguishing A 

We see that w' is related to the nearness of the equa- 
tion of state ratio to the cosmological constant value, i.e. 
1 + w, and the inverse of the characteristic field scale 
of the potential, V^/V (sometimes phrased as a slow 
roll parameter, Mp|V^/y| -C 1). If w is readily distin- 
guished from —1, then we know the dark energy is not a 
cosmological constant, regardless of the value of the time 
variation w' . More difficult is the case where e = 1 + w 
is a small quantity. Then it will be quite important to 
determine whether w' is zero or not. Eq. © guides us in 
following the dynamics in the w-w' phase space. 

The first point to notice is that the reciprocal of the 
characteristic field scale is not generically a small pa- 
rameter useful for a "slow roll" approximation. Figure Q] 
shows curves of constant field scale 



(J) 



in the w-w' plane. Only a tiny sliver of the phase space, 
plus a small hump, satisfies the slow roll approximation; 
unless one is exceedingly close to the cosmological con- 
stant behavior there is substantial dynamics in the field. 

When 1 + w <C 1, the second term in Eq. JfjJ dominates 
and w' > (creating the "humps"), while for less nega- 
tive w the first term can dominate and drive w' < 0. This 
driving occurs closer to w = — 1 for larger Within the 
large-$ hump, the dark energy looks similar to a cosmo- 
logical constant. Suppose one conjectures some physics 
limiting the size of the field scale, equivalently leading to 
avoidance of flatness in the potential. An upper bound 
on <i> in the scalar field behavior would impose a bar- 
rier around the cosmological constant A, saying that the 
scalar field dynamics must be distinguishable from A, if 
it is not A. 

The second panel of Fig. ^ zooms in to illustrate this 
barrier for $ < Mp, which rules out any freezing field, 
and any thawing field with 1 + w < 0.004. So any 
scalar field theory with field scales barred from being 
transPlanckian are distinguishable from A at this preci- 
sion. If the restriction uses, say, <E> < Mp/y/&K (i.e. the 
Planck mass is defined in terms of Newton's constant as 
Gjv = (87rA/p) _1 rather than Gn = Mp 2 as above) then 
the limit becomes 1 + w > 0.1. 

The physical origin for the conjectured limit on the 
characteristic field scale is not clear, but the implications 
are important enough to consider the possibility. Since 
the field scale is related to the inverse of the flatness of the 
potential, then physics that perturbs a flat direction in 
the potential would give this effect. In some supersym- 
metric models, loop corrections generate a logarithmic 
tilt V ~ 4> n hy{<f)/n) HH. This would give $ » <j>/n 
(or <j) ln(^//i) for n = 0), and restricting <p < Mp (for 
the effective field theory to be valid), provides the limit 
$ < O(Mp). However, the generic breakdown of slow 
roll is independent of any $ upper limit conjecture and 
we do not consider the latter further. Hierarchical pa- 
rameters to replace slow roll are discussed in WII 
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FIG. 1: Contours of the inverse of the characteristic field 
scale $ = — V/Vrf are plotted in full and zoomed versions, 
for present matter density f2 m = 0.28. Note that only a slen- 
der crescent of the phase plane obeys the conventional slow 
roll criterion Mp/$ <C 1. The dashed red curve shows the 
null line = 0. Dotted lines show the w ~ — 1 limits of 
the thawing region (upper pair: w' = 3(1 + w), (1 + w)) and 
freezing region (lower pair: here w' — w(l + w),3w(l + w); 
not shown in bottom panel). If physics limits the characteris- 
tic field scale to be less than the Planck mass, then a barrier 
forms (second panel) around the cosmological constant, al- 
lowing only models in the shaded part of the thawing region. 
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B. Driving and dragging 

Returning to the Klein-Gordon equation, we can un- 
derstand behavior in the w — w' phase space by first a 
general and then a specific analysis of the terms. Writing 
Eq. as cf) + 3H<j> = —V^, we can require — > 0, 
reflecting that the field rolls down the potential to large 
(j). Using Eq. (jJJ for 0, we write this condition in terms 



distinct thawing and freezing regions, respectively above 
and below this line. 



> 



-3(1 



''). The boundary defines 



of to, w as w' 
the null line V.^ — discussed further below. Similarly, 
writing the Klein-Gordon equation as + = —3H(j> 
and again using that the field rolls to larger values with 
time implies —3H(j> < 0. This reduces to the condition 
w > -1. 

If we flip the direction of the field motion and potential 
slope, i.e. the field rolls down the potential to smaller 0, 
then these conditions remain. (I.e. V — is still nega- 
tive, and so the transition from Eq. (JSJ to Eq. © flips the 
sign of the second term on the right hand side of Eq. JJjJ , 
canceling the reversed sign of V^.) However, if we make 
the field move up the potential then we have the relations 
w' < —3(1 — w 2 ) and w < — 1, in the phantom region, 
i.e. the boundary lines just continue smoothly through 
w = — 1. 

Finally, we can move to the right hand side to obtain 



3#0 + V A = 



This divides the 



plane into re- 



gions where the field is accelerating or decelerating, with 
the boundary being the coasting behavior = 0. This 
line corresponds to the condition w 1 = 3(1 + w) 2 . Larger 
values of w' arise from a field accelerating down the po- 
tential, while smaller values come from a field deceler- 
ating (this motion of the field should not be confused 
with the accelerating expansion of the universe, which 
can hold for either region) . These three boundaries - the 
null line w' = —3(1 — w 2 ), coasting line w' = 3(1 + w) 2 , 
and phantom line w = —1 - define the broad character- 
istics of the phase plane. 

To investigate the dynamics further, we must examine 
the dominance of the different terms in the equation of 
motion . The first term is the acceleration of the field, 
the second term a friction term, or Hubble drag, due to 
the expansion of the universe, and the third is a driving 
term due to the steepness of the potential. We define 



X 



Y = 



H<j> 



V, 



(8) 
(9) 
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FIG. 2: The physics behind the dynamics becomes more ap- 
parent when looking at the ratios of the terms in the Klein- 
Gordon equation of motion, e.g. field acceleration vs. friction 
vs. potential driving. The key general features of the phase 
plane are the null line = 0, the coasting line = 0, and 
the phantom line w = — 1. Thawing (with w' > 0) and freez- 
ing (with w' < 0) regions are bounded by green, dot-dash 
curve pairs and are defined through the physical dynamics. 



In the accelerating field region, the friction term is 
the major determinant of behavior initially, as the field 
evolves away from a frozen (cosmological constant-like) 
state in the matter dominated epoch. The upper bound- 
ary of the thawing region is given by X = 3/2, where this 
value follows directly from the exponent of the expansion 
history, t ~ a 3 / 2 for matter domination. Thus fields that 
thaw during matter domination begin to move along the 
X = 3/2 line (see discussion, and Figure 2, in Paper 1). 

We can translate any acceleration to friction ratio X 
to a w — w' behavior through 

w' = 2X(l + w)+3(l + w) 2 = (l + w)(3 + 2X + 3w) 
m 2X(l + w), (10) 



Figure [3 shows curves of constant X, Y in the w — w' 
plane. Note that X = Y = corresponds to an epoch of 
coasting in the scalar field dynamics, = 0, as discussed 
above. This is nongeneric, as the field would need to be 
finely tuned to neither accelerate due to the slope of the 
potential nor decelerate due to the Hubble drag, but be 
perfectly balanced. Indeed, the dynamics of scalar fields 
in Paper 1 avoid this region, causing the split into the 



where the last approximation is for 1 + w 1. Note 
that the linear boundaries used in Paper 1 were (good) 
approximations to the general parabolic behavior. The 
value X = gives the coasting line and X = — 3 gives 
the null line. 

Thus, the upper thawing boundary X — 3/2 corre- 
sponds to w' w 3(1 + w) for 1 + w <C 1. If the rolling 
field then enters a region where the potential slope is 
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shallower (as usually happens), then the field will accel- 
erate less and curve toward the <f> — line. Since to- 
day the field cannot have rolled so far that f^c > 0.8, 
the dynamical track remains within the thawing region 
1+w < w' < 3(l + w), i.e. 1/2 < X < 3/2. For potentials 
that steepen as the field rolls down, e.g. PNGB models 
with the field starting near the top of the potential, the 
tracks instead lie above the X = 3/2 or w' = 3(1 + w) 
line. The PNGB potential also steepens more rapidly for 
small symmetry scales /, and indeed, as mentioned in 
Linder ,5j, PNGB models roughly follow w' — F(l + w), 
where F is proportional to the inverse of /. (Also see 
@ for discussion of PNGB models, fine tuning, and slow 
roll.) 

In the decelerating field region, the steepness of the 
potential impacts the freezing. As the potential becomes 
shallower, the friction is more effective. In the limit of a 
flat potential, one obtains the dynamics track given by 
the null curve, = 0, in Fig. [21 As given in Paper 
1, this corresponds to w' = -3(1 - w 2 ), and the skating 
model of 0,0- In terms of the friction and driving terms, 
X = — 3 and Y = oo. The lower boundary of the freezing 
region lies along the Y = 1 (X = —3/2) line, equivalent 
to w 1 = Sw{l + w). (See WI Bl for a rationale.) We will 
later see that this line also has physical significance. 

The general relation between the acceleration to steep- 



ness ratio Y and the 



']' = 3(1+10) 



1-Y 
1 + Y 



track is 



= 3w{l + w) + 3{l+w) 



1 -Y 



1 + Y 
(11) 

Again we have a parabolic behavior. The coasting line 
has Y = 0, and the null line corresponds to Y — > oo. One 
could use either variable X or Y throughout the phase 
plane, since 



X 



Y 



1 + Y 



Y 



X 



X 



(12) 



but this somewhat obscures the physics of friction and 
steepness. (That said, we note that the thawing/freezing 
boundaries are fairly symmetric in X, with the outer 
boundaries at X = ±3/2 and the inner ones at X w 
il/2-) 1 

Note that it is not legitimate to assume tracking be- 
havior (where the equation of state is constant and re- 
lated to the dominant component's equation of state) to 
impose limits on regions in the phase plane, as for exam- 
ple [3, El did to try to tighten the constraints of Paper 
1. For one thing, not all freezing models need start as 
trackers. Secondly, just because the most negative value 
of w' lies above some boundary curve does not ensure 
that the entire trajectory remains above the curve. Most 



1 Recall from Paper 1 that the upper bound on the freezing region 
is not sharply defined, and extends somewhat above the w' = 
w(l + w) line shown for convenience in this paper. 



importantly, the tracking approximation breaks down as 
the dark energy first becomes significant, so it is inappli- 
cable for much of the observable dynamical history. 



III. MORE SPECIFIC DYNAMICS 

In addition to analyzing general behavior through the 
Klein-Gordon equation terms, we can investigate the 
properties of the phase plane or specific track families 
in terms of other variables. These could include work- 
ing from the equation of state w — w' relation directly 
or the cosmic expansion history a(t). While not quite 
as insightful as the physics motivated driving and drag 
terms, they can highlight interesting properties. 



A. Mocker models 

Consider a model with dynamics given by w' — Cw{l + 
w). Note that while such an equation forms the bound- 
aries of the freezing region, freezing models do not follow 
such a trajectory but rather are almost orthogonal to 
such tracks (at least initially). So we are talking about 
fundamentally different models. The behavior of the dark 
energy equation of state and density follow 



1 __» L _ a a 



w(a) = -1 + 

l + w 

Pde(a) = Pde(l) [(1 + Wq) 



C 1 3 / C 

a — wq\ 



(13) 
(14) 



where wq is the equation of state today. 

In the past, a <C 1, the component will act like addi- 
tional nonrelativistic matter, with w = 0, while in the 
future it will approach a cosmological constant. Since 
such dark energy sometimes looks like dark matter and 
sometimes like a cosmological constant, we call this a 
"mocker" model. These are basically what are known as 
"quartessence" models (see for an overview) , of which 
the Chaplygin gas is one example. However, we develop 
them directly from the phase space dynamics rather than 
an ansatz for the pressure, so the dynamical behavior is 
more general. For example, a constant pressure model 
could be a cosmological constant, or could be a mocker 
with C = 3. Note C = 3 gives precisely the expression for 
the lower boundary of the freezing region (or is Y = 1 
or X — —3/2 in the notation of £111 Bp . We name the 
w' = 3w(l + w) curve the constant pressure line (also see 
[Io| in the context of barotropic fluids) . 

Such combined behavior models are heir to all the 
usual problems of trying to unify dark matter and dark 
energy, e.g. growth instabilities of density perturbations 
|12l Il3j . Analysis of perturbations requires knowledge 
of the full theory, however. Merely from the phase plane 
dynamics, though, we can see trouble arising for such uni- 
fied models. As C gets smaller, the model moves along 
its trajectory more quickly, acting less like dark matter 
except at very early times, 1 + z ^> [— Wo/{l + wq)] 1 / . 
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Conversely, as C gets larger, acceleration of the cosmic 
expansion occurs later and the model becomes a poorer 
fit to a host of (purely geometric) observations such as 
supernova distances and the distance to the CMB last 
scattering surface. 



B. Relation to parameterized w(a) 

The approach taken in this article is to examine dark 
energy dynamics directly in the phase plane w—w', where 
a time variable runs along each trajectory. It is useful to 
see the relation of standard parametrizations in terms of 
the temporal behavior, i.e. w(a), to this approach. 

The standard two parameter function w(a) = Wq + 
w a (l — a) was shown by Linder [Til ITsT ] to provide an 
excellent approximation to exact solutions of the Klein- 
Gordon equation in a wide variety of models. In this 
ansatz we have w' = —aw a = w — w^, where the high 
redshift equation of state = wo + w a . This describes 
a straight line of slope 1 in the w — w' plane, and can be 
rewritten as as w' = (1 + w) — (1 + w^). In particular, 
if Woo = — 1 we have exactly the behavior of thawing 
models (lying along lower bound of that region). From 
Fig. 2 of Paper 1, we see as well that many tracking 
models that fit present data (i.e. f^c ~ 0.7 and w < 
—0.8) are reasonably well described, on average, by a line 
of slope unity. Of course this approximation will break 
down in the future, as the field freezes more fully, turning 
toward the cosmological constant; at the same time this 
w(a) ansatz loses validity as it moves toward ever more 
negative w. However, since data only exist toward the 
past, we see why the w a parametrization is an excellent 
approximation. 

To retain boundedness for both the past and present, 
as well as to allow more dramatic dynamics (essentially 
slopes other than unity inw-to'), one could use the "e- 
fold" model of 16] or "kink" model of H3- Both utilize 
four parameters for their description. The e-fold model 
has a more transparent translation to w' = dw/dlna 
since it also uses dynamics in terms of In a: 



w(a) = Wf + 



Aw 



1 + (a/a t y ' 



(15) 



where r is the transition rapidity in units of units of 
e- folds In a, a t is the transition scale factor, Wf is the 
asymptotic future value, and Aw = w p — Wf is the dif- 
ference between asymptotic past and future values. 
In the w — w' phase plane we have 



w 



-t(w — Wf) I 1 



Wf 



Aw 



(16) 



Note that as we found in the Klein-Gordon equation 
analysis of £111 Bl the equation for w' is quadratic in w. 
We can identify several special cases. If the asymp- 
totic future state is deSitter (wf — —1), then w' = 
t(Aw _1 -1)(1+w) + (t/Aw) io(i+«;). This looks like the 



sum of a thawing model and a model in the freezing re- 
gion, i.e. the dark energy can be viewed as the sum of two 
components. If we further take Aw = 1, then we remove 
the thawing component and end up with w' — tw(1 + w) 

- a mocker model with w p = and Wf = — 1. 
Starting instead with an asymptotic past state of w p — 

— 1 gives w' = t (1 + w)(w — Wf)/(— 1 — Wf). In the 
limit Wf — > oo (i.e. not worrying about the region where 
there is no data) this gives a thawing model w' = r(l + 
w). Thus the four parameter e-fold ansatz is also quite 
versatile. The rapidity parameter is directly related to 
both the slope of and the velocity along the phase space 
trajectory, and ties in with the steepness of the scalar 
field potential, as we saw in i jllBl with the PNGB models 
where the slope was proportional to the inverse of the 
symmetry scale /. 

Finally, one could invert the situation and go from a 
parametrization in the phase plane to derive the function 
w(a). For example, a track w' — A(l + w) + B(\ + w) 2 , 
which we have seen is a common form, implies 



l+io= (1 + w )/ [(l + a;> 



—A 



(17) 



where x = (B/A)(l + wq) defines the present along the 
trajectory (equivalently the dimensionless matter density 
Q, m today). Note that while the trajectory has two pa- 
rameters, the equation of state has three parameters since 
we must have a parameter running along the track. At 
high redshift, if A > then w — * — 1 and we have a thaw- 
ing model, asymptotically independent of B. If A < 
then w(z > 1) = -1 + (-A/B) and w' — > 0, i.e. it be- 
gins like a tracking model. It reaches a minimum w' min = 
-A 2 /(4B) at w if = -l- A/(2B) = [-1 + w(z > l)]/2; 
that is, the trajectory is a parabola from its tracking 
value of the equation of state to its future, cosmologi- 
cal constant value. A mocker model is the special case 
C = —A = B. For completeness, we give the dark energy 
density: 



Pdc(a) = Pde(l) (1 + x - xa ) 



A\Z/B 



(18) 



Of particular interest is the "leveling" model where, in 
loose physical analogy to the inflationary power spectrum 
tilt n — 1 being driven to zero by large numbers of e- 
folds of expansion, the equation of state tilt 1 + w is 
driven to zero by the deSitter expansion as the energy 
density approaches a certain constant value, pf. That is, 
take 1 + w = D[pdc(a) — pf]. This is equivalent to the 
above parabolic model with A = —3Dpf and B = —3. 
Another interesting parabolic track is the coasting line 
w' = 3(1 + w) 2 . This corresponds to not a leveling but 
a tilting, with 1 + w = (1 + wo) [pde(l) / Pde(a)] , so w is 
tilted away from —1 as the energy density decreases. 



C. Acceleration and jerk 

One could leave behind the physics of the accelerating 
phenomenon and instead use variables in terms of the 
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acceleration itself, though this seems less appealing. The 
deceleration parameter 



-ad/a 2 



1 



-wCl de (a) 



(19) 



and the jerk 



j = a 2 a/a 3 = 1 - -Vl de {a) [w' - 3w(l + to)]. (20) 

We also have j = q+2q 2 —q'. Note that interpreting q and 
j or q' as Taylor expansions of the expansion is of strictly 
limited use (since observations span A In a ~ O(l)) and 
can be dangerous 0] . 

Furthermore, there is the same ambiguity there was 
with using pressure as a variable. We note that any 
model where it touches the constant pressure line w' = 
3w(l + w) has j — 1; this is equivalent to X — —3/2 in 
the notation of £111 Bl The two standard special cases of 
j = 1 lie at the ends of this line: an Einstein-de Sitter 
universe with w — "dark energy" and a ACDM universe 
with cosmological constant dark energy. 

The constant pressure line is also related to the adia- 
batic sound speed of the dark energy. (Note this is not 
the true sound speed of perturbations arising from the 
microphysics of whatever the dark energy is.) The adia- 
batic sound speed 



V 1 
- = w — — 



3 1 + w ' 



(21) 



and we see it vanishes for w' — 3w(l + w). On the null 
line, the adiabatic sound speed equals the speed of light 
(the same as the true sound speed for a canonical scalar 
field). Models below the null line would need to have 
cl > 1. 



IV. BEYOND SCALAR FIELDS 

For the cosmic expansion dynamics we can always de- 
fine an effective equation of state even if the accelerating 
mechanism is not a scalar field |l|j, through 



1 din 6 H 2 

WeS = -1 - 7; ' 

6 dm a 



(22) 



where 6H 2 = (H/H ) 



,-3 



is the unknown part of 



the Hubble parameter, that not due to matter. So it is of 
interest to see to what extent the dynamical behaviors we 
have discussed carry over to the w e g — w' eS plane. That 
is, are freezing and thawing behaviors more general than 
for scalar fields, and do the null and coasting lines still 
play a role? 

Due to the diversity of possible accelerating physics 
we do not present a general analysis of these important 
questions but rather calculate some specific cases. 



A. Scalar-tensor gravity 

Scalar-tensor theories modify the Einstein- Hilbert ac- 
tion with both an additional scalar field and a coupling to 
Ricci scalar curvature R. These are of great interest as a 
comparison in tests of general relativity and also because 
gravitational theories involving a nontrivial function of 
the scalar curvature can be transformed to scalar-tensor 
theories. See for a general introduction. 

We consider coupling of a general form in the scalar 
field, but linear in the curvature. So the general rela- 
tivistic R/(8ttG) — > F{(f))R. One then obtains the usual 
Friedmann expansion equations, with extra terms giving 
an effective dark energy density [2(j 

Pst = V(4 ) ) + (l/2)H 2 (q-l)(q+5)F 2 +3H 2 [(87TG)- 1 -Fl 

(23) 

where q is the deceleration parameter, F$ = dF/d(j), and 
the last term involves the change of the gravitational 
strength from Newton's constant G. 

From this density one can then define 8H 2 = 
8itGpst/(3H 2 ) (note that we use the usual G here since 
the deviation is absorbed into pst as mentioned above). 
From this the equation of state w e g and its variation w' eS 
can be calculated. A key quantity will be F/F£ = ujjbd- 
This is the Jordan-Brans-Dicke parameter and its inverse 
must be very small according to solar system tests. Ex- 
panding F about the present, 

F(a) fa (87TG)- 1 - (1 - a)F^/{aH) (24) 



(8^G)- 1 -z(l-g)^ 2 . 



(25) 



So the ratio of the second to the first term (~ ^jb D ) is 
small, and gravity is nearly Einsteinian. 

But this means that the first term in pst dominates 
(unless 8nGV/H 2 <C 1, but then it doesn't affect the 
expansion and there is no acceleration). Thus, the re- 
striction of the scalar-tensor theory by solar system con- 
straints means that its effective equation of state must be 
very close to a cosmological constant - within ~ ^jbd- 
Since ijJj BD < 2.5 x 10~ 5 , this would be rather chal- 
lenging to distinguish from a cosmological constant with 
cosmological observations! One possible loophole is if the 
solar system limits on the scalar coupling should not be 
applied to a cosmological situation because of the dif- 
ferent spacetime backgrounds with very different scalar 
curvatures. This arises for example in chameleon sce- 
narios 21] . The most stringent cosmological bounds on 
varying G arise fro m p rimordial nucleosynthesis and give 

^JBD ~ 3 x 10 " 3 ( but see E3)- 

Calculation of the effective phase plane parameters 
finds 



w cS (z 

w 'eff( z 



0) 
0) 



-1 + 0A6/LJJBD 
-0.36/lujbd- 



(26) 
(27) 



While even with only the cosmological bounds on ujjbd 
these are quite close to the cosmological constant in the 
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phase plane, it is interesting to note that the current 
values lie along w' — 0.78w(l + w), in the freezing re- 
gion. Indeed its trajectory is a freezing one, with scalar- 
tensor theories asymptotically attracted to general rela- 
tivity |26|. and to w = — 1. One last thing to note, 
however, is that because scalar-tensor theories possess 
anisotropic stress, the growth of density perturbations 
will be modified from the quintessence case (see, e.g., 
[21} and references therein). 



B. Braneworld cosmology and H c 



In a braneworld cosmology 28, 29], effective accelera- 
tion appears due to a weakening of gravity on large scales 
as it "leaks" from our brane into a higher dimensional 
bulk. The Friedmann expansion equation becomes 



H 2 - H/r c = 87rGp m /3, 



(28) 



where r c is the crossover distance and p m the matter den- 
sity. The effective equation of state due to the modifica- 
tion is w e g = —[1 + n m (a)] _1 [3(j. Its trajectory in the 
phase plane is plotted in Fig. El Note that it looks like a 
freezing model, and will indeed approach a cosmological 
constant in the asymptotic future. 

The position along the trajectory is a time variable, so 
taking the present to be, say, Jl m = 0.2 would extend the 
solid curve in Fig. [21 slightly further (since the figure uses 
J1 TO = 0.3). We also see why it is so well approximated by 
a,wo~w a model, as discussed in fllll Bl (Recall, however, 
that w a is actually defined at z = 1, not z = 0, to give 
the best physics fit [14(.) 

We can further generalize the modification to SH 2 ~ 
H a H3- Then we find 



a 



-Jim (a] 



v' cS = 3w(l +w)[l -(2/a)(l +«;)]. 



(29) 
(30) 



Recall that the braneworld case above corresponds to 
a = 1. For acceleration today (with Jl m = 0.3), we 
require a < 1.57; for w < —0.8 today we require a < 
0.91. Note that all H a modified gravity models will look 
similar (one does require a < 2 for a negative equation 
of state at early times). See WIBI for discussion of their 
tracking behavior. In particular, they all approach the 
cosmological constant along w' = 3w(l + w), what was 
called the constant pressure line for scalar fields. Their 
tracks must always lie between w' = 3u>(l + w) and the 
w' = axis. When a < 0, the trajectories switch to the 
phantom regime with w < — 1, but the bounds still hold. 



V. POLYTROPIC DARK ENERGY 

An interesting, if phenomenological, way of obtain- 
ing acceleration is to modify the Friedmann expansion 
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FIG. 3: The braneworld model of modified gravity behaves 
like a freezing scalar field. The solid black line shows its 
trajectory: starting at high redshift at w — —1/2, w' — 0, 
with the red crosses indicating z = 3, 2, 1 before the line ends 
at the present. The dashed blue line shows its future behavior, 
approaching a cosmological constant, and lying within the 
quintessence freezing bounds of w' = 3ui(l + w) and w' = 
w(l + w) given by the dotted green curves. A model with 
a generalized, extra H a term (for the braneworld, a — 1) in 
the Friedmann expansion equation looks similar; the magenta, 
dot-dash curve is for a = 0.5. 



equation but keeping a pure matter universe. While this 
leaves open important questions about its relation to fun- 
damental theory and the growth of density perturbations, 
we can investigate some general aspects of the effective 
equation of state dynamics. 

Consider general functions of the matter density 
(sometimes referred to as barotropic models [lO^ 



H 2 = (8ttG/3) g(p) = (8ttG/3) [p + f(p)}. 



(31) 



The quantity f(p) will act like an effective dark energy. 
Using Eq. (1221) we can define 



WcS 



1 



din/ 
dhip 
d 2 In/ 
d In p 2 ' 



(32) 
(33) 



and identical relations hold for the total equation of state 
u>tot and its variation w' tot upon substituting g for /. 

The simplest example of such a modification is / ~ p n , 
the Cardassian model of [3^1 • From the above equations 
we see that it corresponds to a constant equation of state 
w = — 1 + n. If we require w < —0.9 (as observations 
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favor for a constant equation of state), then n < 0.1; 
unfortunately p 1 / 10 does not obviously appear to be a 
natural modification of the Friedmann equation. 

We can investigate the phase space dynamics by relat- 
ing w' to w. 

o 2 d 2 f 

w' = 3 W (l + W )-3^54. (34) 
/ dp 2 

Immediately we see that whether the effective dark en- 
ergy lies below the freezing region or not depends on the 
sign of d 2 f /dp 2 . An equivalent question is whether the 
effective pressure is decreasing or increasing with time 
(since the w' = 3w(l + w) line is that of constant pres- 
sure) . 

The model will follow the mocker model track w' = 
3u>(l + w) if / = A + Bp. This form is equivalent to a 
redefinition of Q m , e.g. f2 m — > fi m (l + B), plus a cosmo- 
logical constant A. As such it has a nonzero minimum 
in its effective potential, allowing it to reach the freezing 
boundary. 

It is important to remember that the analysis in this 
paper and Paper 1 applies to the dynamics of the dark en- 
ergy itself. Trajectories in the iwtot — w' tot plane convolve 
the matter and dark energy components, mixing the dy- 
namics and so not giving rise to the clear differentiations 
and regions found. This is why Ref. [10| appears to find 
a violation of the freezing bound and even null bound 
for some barotropic models; those are actually phantom 
models in the dark energy phase space, but are dragged 
by the matter to lOtot > — 1- We make this more explicit 
later in this section. 

For a richer dynamical behavior we propose a class of 
modifications of the Friedmann expansion equation we 
call polytropic models [33]. Here 

H 2 /(8wG/3) = g(p) =p[l+ (p/p^f. (35) 

At densities much greater than some crossover value p*, 
e.g. at high redshift, the Friedmann equation is standard. 
At low densities, the expansion is modified, with w to t 
asymptotically approaching —af3. If we want a future 
deSitter state, we could choose (3 = l/a. For just the 
effective dark energy equation of state, the value in the 
past is w c s — —a, and in the future of course it dominates 
so w e a = —a/3. 

Figure 0] illustrates the phase plane dynamics. The 
first panel takes (3 = l/a, so that the future state is de- 
Sitter. Models with a < 1 have w > — 1 and act like 
freezing models, starting from a constant w = —a and 
today (marked by crosses) lying in the freezing region, 
before heading toward the cosmological constant. Phan- 
tom models have a > 1 and act like mirror images of 
freezing models, even to lying within the phantom freez- 
ing region today. 

In the second panel, we fix f3 — 1/2. While the models 
start at the same phase space point as the previous mod- 
els with the same a, now their endpoints are different. 
Indeed for a < 2/3 the acceleration of the expansion is 
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FIG. 4: Polytropic models, with /3 = l/a in the first panel 
and with j3 = 1/2 in the second panel. Curves are labeled 
with a and crosses show the present. Dotted curves show the 
canonical scalar field freezing region (first panel) or thawing 
region (second panel). In the first panel, models evolve toward 
a cosmological constant, within the freezing region. In the 
second panel, the asymptotic state is not deSitter (except for 
the a = 2 curve, which is the same in both panels) and the 
dynamics is distinct from a canonical scalar field. 



a temporary phase. Furthermore, the trajectories with 
a < 1.5 do not lie in the freezing regime and both reg- 
ular and phantom models have w' > 0. Such polytropic 
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models without a deSitter future will be clearly distin- 
guishable from both freezing and thawing quintessence. 

Can we give some physics motivation for the polytropic 
form, aside from its simplicity and proper asymptotic be- 
havior? When (3 — 0, there is no modification; when 
(3 = 1 we have the power law modification of the Car- 
dassian case, with n = 1 — a, hence a constant w = —a. 
There are some motivations for power law modification 
from higher dimension theories (for n < 1 see j^l, the 
nonaccelerating n = 2 arises in Randall-Sundrum brane 
scenarios [3_^). When (3 = 1/2, the modification is sim- 
ilar to that from a Chaplygin gas (35|, as we see below; 
this has claimed motivation from Born-Infcld actions and 
brane solutions [3(|. So at least the polytropic form uni- 
fies different prescriptions for dark energy. 

Note that as (3 increases from zero, for fixed a, the size 
of the "hump" in the trajectory decreases and the future 
value of w moves back toward the initial value. At (3 = 1 
the trajectory collapses to a point at w = —a, w' = 0. 
For even larger /3, the hump is flipped (i.e. the sign of 
w' changes) and again increases in size, with the future 
value of w drawing away to more negative values. 

When we plot the same models as in the first panel 
of Fig. but in the ui to t — w' tot plane, in Fig. [S] we see 
that the models that were phantom in the effective dark 
energy lie in the region w' tot < 3u>tot(l + wtot)- Further- 
more, when a > 2 they can even lie below what was 
the null line w' = —3(1 — w 2 ). In a nice analysis of 
barotropic models, Scherrer 10] noted something simi- 
lar (his barotropic models are a function of an arbitrary 
component density not necessarily matter density). For a 
barotropic perfect fluid the adiabatic sound speed is the 
physically relevant sound speed for perturbations (but 
not in the quintessence case, or in a general multicompo- 
nent case), and Scherrer's bound of w' < 3w(l + w) holds 
for c 2 > (cf. Eq. [71] here). In general, however, this is 
not a violation of the bounds of this article and Paper 
1 because it occurs only when the adiabatic assumption 
holds, e.g. when viewing the total equation of state, not 
the properties of a (non-adiabatic) effective dark energy. 

Suppose, however, we chose to fit the dark energy com- 
ponent itself, rather than full energy density entering the 
Friedmann equation, by the polytropic form Eq. I|35|) . 
This is somewhat strange to do, since then the effective 
dark energy contains a matter-like part, in addition to 
the pure matter density, and the polytrope was designed 
to be the modified Friedmann equation as a whole. If 
we do so, though, then the dark energy equation of state 
phase plane (and not the total equation of state) is rep- 
resented by the curves in Fig. [3] Moreover, the curve 
with a = 2, /3 = 1/2 is the trajectory of the Chaply- 
gin gas. The generalized Chaplygin gas with pressure 
p ~ — p~ a s"3 corresponds to polytropic dark energy with 
a = a gcg + 1, P = l/a- Whenever a(3 = 1 (cf. |37jl 
we have mocker behavior with w = 3aw(l + w). As 
stated above, however, taking the dark energy itself to 
be polytropic means hiding both a cosmological constant 
(if (3 = l/a) and a spurious matter density within the 
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dark sector. 



VI. HIERARCHY PARAMETERS 

Returning to more general dynamics, we saw in fJH] 
and in particular from Fig. ^ that we do not have the 
standard inflation slow roll perturbative expansion pa- 
rameter. Explicitly, we do not have small V^/V in the 
freezing or thawing regions unless w % —0.995; more gen- 
erally, —MpV^/V > 5\/l + w today. We examine here 
whether we can substitute a physics based hierarchy of 
dynamics parameters. 

A. Slope parameters 

The Klein-Gordon equation can be rewritten dimen- 
sionlessly as 

ft' + (2-q)<f)' = -V^/H 2 = m , (36) 

where as before a prime denotes a derivative with respect 
to In a, and q = — ad/ a 2 is the deceleration parameter. 
In terms of the w — w' dynamics equation, 

w' = -3(1 - w 2 ) - V2(l - w 2 ) V^/{HV 1/2 ) (37) 
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-3(1 -w 2 ) + ^2(1- w 2 )r ]2 . 



(38) 



So the parameters rji, ij 2 are called out by the physics. 
If they are small, one could solve the equations pertur- 
batively. Note that 772 = r?i(V '/ 'H 2 )^ 1 ^ 2 so we always 
have 772 > T]i . We violate the lower bound of the freezing 
region, w' < 3w(l + w), when 



V2 < 



9 



22 



where e = 1 + w is the tilt parameter. The analogous 
condition such that w' < is rj 2 < (9/2)e(2 — e), and 
such that field is decelerating (4> < 0) is rfe < 18e/(2 — e). 

Neither 771 nor 772 are particularly small unless e = l+w 
is. For example 772 > l.5y/e, giving 772 > 0.1 for w > 
-0.995. Even 771 > 0.1 for w > -0.94. The hierarchy 
among the parameters is fixed for the region of interest: 
e < 771 < 772, so there is no phase space classification in 
this respect, as there is for large field, small field, and 
hybrid models in inflation [3^]. However that hierarchy 
involved the second derivative of the potential, so it is 
worth a brief look at that quantity. 



B. Tracking parameter 

The tracking parameter is defined to be 



vt 



(40) 



This is not generally a small parameter cither. Indeed, 
models whose energy density tracks |39j the evolution of 
the dominant energy component fulfill the conditions 



r > 1 



dln(r- 1) 
din a 



< 1. 



(41) 



Within the class of tracking models (so now a particu- 
lar subset of scalar field cosmologies), at high redshifts 
within the matter dominated epoch the field obeys 



1 



w 



2(1 + w) 



(42) 



(Note the equation of state deviation 1 + w at high red- 
shift may not be small.) 

In 9IVBI we found that H a models (including 
braneworlds) follow freezing trajectories. This is not sur- 
prising because they are basically trackers. A component 
starting with constant w at early times is equivalent there 
to a modification H a with a = 2(1 + w). The tracking 
parameter T — (2+a)/(2a) and it initially acts like an in- 
verse power law potential V ~ <\>~ n with n = 2a/ (2 — a). 

To relate the dynamics of the time variation to the 
tracking condition, we invert Eq. (|42|l to write 



w 



w 



-2(r- 

dw 
din a 



= w(l - 



2(r-i)] 

dln(r- 



w 



G?lna 



(43) 
(44) 



However, this is of limited use since we are unlikely to 
be able to probe w' in the high redshift, z » 1, regime 
where tracking might hold. Strong acceleration today, 



with 



< 



-0.7, requires the breakdown of tracking. 



However, the analogy to H a models presents an im- 
portant insight into why the general freezing region is 
bounded below by w' = 3u>(l + w). At early times the 
contribution to the Friedmann expansion equation by the 



(39) dark energy is small and SH 



-3(l+«<) ^ ff2(l+v 



That is, the effective a ~ 2(1 + w). One can generalize 
this to a = 2(1 + w) when time variation of the equation 
of state becomes relevant, where angle brackets denote 
an averaging over In a. At late times, the dark energy 
density dominates the expansion, 5H 2 ~ H 2 , as it ap- 
proaches w = — 1 (freezes). For any epoch we can define 
an instantaneous value of a. Equation (|30|l then gives 
the relation for w' . As freezing models approach w = — 1, 
Eq. (|30|l indicates they should do so along w' — 3w(l+w). 
Furthermore, since the bracketed term is less than one, 
then this trajectory represents a general lower bound to 
the freezing region. 2 



C. Dynamics, Mass, and Spatial Inhomogeneities 

Dynamical models must also possess spatial inhomo- 
geneities in the field at some level. The equation for 
these is given by perturbation of the Klein-Gordon equa- 
tion 0, 



■ + 3H5(j) + {k 2 + V^Sc/) = -h(f>/2, 



(45) 



where k is the wavenumber and h is the trace of the 
metric perturbation [40]. Just as matter density per- 
turbations are damped on scales below the Jeans length 
related to the sound speed in the background medium, 
so the spatial inhomogeneities in the scalar field will be 
absent on length scales less than that corresponding to 
the effective mass WV^. 

Using eq. © and = V<j,/4> we can calculate the 
critical mass scale (also see [ill]), with 



V^/H 2 = (2 + 3w + 2q) 



1 



1 

4 V 1 



1 w" 3 
"2lT^ + 4 (1 " W)(5 + SW + 2q)m 

Note w'/(l + w) and w" /(\ + w) are well behaved and 
generally nonzero as w — > — 1. 

This shows that the goal of exploring the temporal and 
spatial dynamics of dark energy runs into double jeop- 
ardy. If the time variation is weak, \w'/(l + w)\ <C 1, 



2 One caveat involves dark energy models that possess an internal 
cosmological constant, i.e. nonzero minimum to the potential, 
or otherwise act as the sum of two or more components. These 
cannot be represented as H a models and the freezing bound does 
not apply. 
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then the effective mass m = y/V^ < H. (The scale H 
today corresponds to 10 -33 eV; dark energy would be a 
very light scalar field). This means the Compton wave- 
length of the scalar field perturbations is larger than the 
horizon, and so spatial inhomogeneities are also difficult 
to detect. Thus, for \w'\ < 1 + w, i.e. between the thaw- 
ing and freezing regions there is a "dead zone" of phase 
space, where we can detect neither time variation nor 
spatial inhomogeneity. 

For appreciable time variation, \w'\ > 1 one can 
have m > H and so the possibility of subhorizon cluster- 
ing. However for models within the freezing or thawing 
regions, one is restricted to m < 2H so this could only 
occur on the largest scales (largest angles or lowest mul- 
tipoles). Note that as the field approaches w = — 1, the 
mass stays nonzero (except it vanishes along the upper 
boundary of thawing and along the null line). However, 
the amplitude of spatial perturbations vanishes as can 
be seen from Eq. (|45|l with <j> = 0. These properties 
make scalar field inhomogeneity an extremely begrudg- 
ing probe of the nature of dark energy, much less friendly 
than the dynamics. 

VII. CONCLUSION 

We have deepened and elaborated the understanding 
of the role that the dark energy dynamics, through the 
w — w' phase plane, can play in leading our understand- 
ing of the nature of dark energy. This includes the foun- 
dations of the null line, coasting line, constant pressure 
line, and phantom line dividing the phase plane into dis- 
tinct, physical regions. We also elucidate the uppermost 
and lowermost boundaries of the thawing and freezing 
regions. 

The physical structure has been extended beyond 
canonical scalar fields, including specific instances 
of modified gravity scenarios such as scalar-tensor, 
braneworld, and H a models, and barotropic and poly- 
tropic generalizations of the Friedmann equation. We 
outlined similarities and differences with the scalar field 
case, showing that many act as freezing fields, and that 
we should be able to clearly distinguish certain models 
that do not possess a deSitter future. Mocker models, im- 



plementing a unification of dark matter and dark energy, 
were shown to have difficulties purely from dynamical 
considerations, in addition to their problems in structure 
formation. 

Dark energy is demonstrated to be generically not 
amenable to a slow roll description - a major difference 
from early universe inflation - as one of its "Goldilocks" 
conundra. This makes the dark energy problem in some 
sense even more challenging than the early universe. 
However, it also opens the possibility that if some phys- 
ical bound can be placed on the flatness of the poten- 
tial, e.g. due to quantum corrections, then this implies a 
barrier around the cosmological constant A model. This 
would offer hope, possibly accessible to next generation 
experiments, that dark energy could definitely be distin- 
guished from A, if it is not A. That would be exciting! 

The dynamics of the dark energy, in the form of the 
equation of state ratio w and its time variation w', pro- 
vides powerful insight into the new physics behind cosmic 
acceleration. Spatial inhomogeneities in the dark energy 
are seen to be much weaker and less forthcoming, unless 
one entertains direct couplings. While we are not guar- 
anteed to zero in on the physics - there is a dead zone of 
minimal dynamics and possibly a "confusion" zone near 
the cosmological constant any highly precise and accu- 
rate result would be an enormous success in enlightening 
us on the dark universe. 
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